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annihilation at 

high luminosity factories using the radiative return method is motivated 
and discussed. A Monte Carlo generator which simulates the radiative 
process e + e~ — > 7 + hadrons at the next-to-leading order accuracy is pre- 
sented. The analysis is then extended to the description of events with 
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1. Motivation 

Electroweak precision measurements in present particle physics provide 
a basic issue for the consistency tests of the Standard Model (SM) or its ex- 
tensions. New phenomena physics can affect low energy processes through 
quantum fluctuations (loop corrections). Deviations from the SM predic- 
tions can therefore supply indirect information about new undiscovered par- 
ticles or interactions. 

The recent measurement of the muon anomalous magnetic moment = 
{g — 2)^/2 at BNL [1] reported a new world average showing a discrepancy 
at the 2.6 standard deviation level with respect to the theoretical SM eval- 
uation of the same quantity which has been taken as an indication of new 
physics. For the correct interpretation of experimental data the appropriate 
inclusion of higher order effects as well as a very precise knowledge of the 
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SM input parameters is required. The BNL experiment plans a new mea- 
surement with an accuracy three times smaller which will challenge even 
more the theoretical predictions. 

One of the main ingredients of the theoretical prediction for the muon 
anomalous magnetic moment is the hadronic vacuum polarization contribu- 
tion which moreover is responsible for the bulk of the theoretical error. This 
quantity plays also an important role in the evolution of the electromagnetic 
coupling olqed from the low energy Thompson limit to high energies. In 
both cases the precise knowledge of the ratio 

„. . a(e + e + — > hadrons) , , 

cr(e + e — ► fj, + /j, ) 

over a wide range of energies is required. The measurement of the hadronic 
cross section in e + e~ annihilation to an accuracy better than 1% in the 
energy range below 2 GeV is necessary to improve the accuracy of the 
present predictions for the muon anomalous magnetic moment and the QED 
coupling. 

In this paper, the radiative return method is motivated and described. 
It has the advantage against the conventional energy scan [2], that the 
systematics of the measurement (e.g. normalization, beam energy) have to 
be taken into account only once but not for each individual energy point 
independently. An improved Monte Carlo generator which simulates the 
radiative process e + e~ — > j+hadrons at the Next-to-Leading Order (NLO) 
accuracy is presented. Finally, the analysis is extended to the description 
of events with hard photons radiated at very small angle. 



2. The muon anomalous magnetic moment and the QED 

coupling 

The Standard Model prediction for the muon anomalous magnetic mo- 
ment consist of three contributions 1 

a, = < ED + af d + a» eak . (2) 

The QED correction, known (estimated) up to five loops [4], gives the main 
contribution (see Table 1 for numerical values 2 ). The weak contribution 
is currently known up to two loops [5]. Finally, the hadronic contribution 
consist of three different terms where the leading one comes from the vacuum 
polarization diagram in figure 1 which furthermore is responsible for the bulk 
of the theoretical uncertainty. 

1 For recent reviews see [3]. 

2 N.A. see [27] for a recent evaluation of the hadronic light- by- light contribution. 



3 



a M x 10 



rr 



QED 116584706 ± 3 

hadronic (vacuum polarization) 6924 ± 62 [8] 

hadronic (light - by - light) -85 ± 25 [12] 

hadronic (other higher orders) -101 ± 6 [13] 

weak 152 ± 4 



Table 1. QED, hadronic and weak contributions to the muon anomalous magnetic 
moment. 




vac. pol. light-by-light 



Fig. 1. Vacuum polarization and light-by- light hadronic contributions to the muon 
anomalous magnetic moment. 

While the QED and the weak contributions are well described pertur- 
batively, the hadronic contribution cannot be completely calculated from 
perturbative QCD. The hadronic vacuum polarization contribution [6, 7, 8, 
9, 10, 11] can however be evaluated from the following dispersion relation 

< d (vac.pol.) = f^f r % K(s) R(s) , (3) 

where the kernel K{s) is a known smooth bounded function and R(s) is the 
hadronic ratio which has to be extracted from experimental data. Under 
the assumption of conserved vector current and isospin symmetry r decays 
can be included also for the evaluation of the dispersion integral. Due 
to the 1/s 2 dependence of the dispersion integral, the low energy region 
contributes dominantly to this integral: 70% of the result comes from the 
7T7T channel and, even more, 90% of it is given by the region below 1.8 GeV. 
The precise measurement of the hadronic cross section in e + e~ annihilation, 
specially at low energies, has therefore a crucial importance for the accurate 
determination of the muon anomalous magnetic moment. 

The running of the QED fine structure constant from the Thompson 



4 



Aa(M z ) x 10 4 


leptonic 


314.98 


top quark 


-0.70 ± 0.05 


hadronic 


276.1 ± 3.6 



Table 2. Leptonic, top quark and hadronic contributions to the running of the 
QED fine structure constant at Mz- 



limit to high energies is given by 

a ( s ) = -, n ) A « = Aa lep + Aa top + Aa had . (4) 

1 — /\a(s) 

Again, the Standard Model prediction consists of several contributions. 
Among them, the hadronic one which can be evaluated from the disper- 
sion relation 

aS r, f°° ds ' R ( S ') 

Aa had (s) = -—Re / — y~ . 5 

3vr J im 2 s' s s if] 

The leptonic [14] and the top quark [15] contributions are well described 
perturbatively while the hadronic contribution [16] gives the main theoreti- 
cal error (see Table 2 for numerical values). The dispersion integral grows in 
this case only as \/s which means that also the high energy points become 
relevant. The region below 1.8 GeV contributes only to 20% of the integral. 
A better knowledge of the hadronic cross section in this region would be 
nevertheless also important to reduce the error. 



3. The radiative return method and the hadronic cross section 

The radiative process e + e~ — ► ^y+hadrons, where the photon is radiated 
from the initial particles (initial state radiation, ISR), see figure 2, can be 
used to measure the hadronic cross section <r(e + e - — > hadrons) at high 
luminosity electron-positron storage rings, like the ^-factory DAPHNE or at 
.B-factories, over a wide range of energies. The radiated photon reduces the 
effective energy of the collision and thus the invariant mass of the hadronic 
system. This possibility has been proposed and studied in detail in [17] (See 
also [18]). A Monte Carlo generator called EVA [17] which simulates the 
process e + e~ — ► 7r + 7r~7 was built. The four pion channel was considered 
in [19]. 

Radiation of photons from the hadronic system (final state radiation, 
FSR) should be considered as the background of the measurement. One of 
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ISR FSR 



Fig. 2. Initial state radiation (ISR) and final state radiation (FSR) in the annihi- 
lation process e + e~ — > 7 + hadrons. 
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Q 2 (GeV 2 ) 

Fig. 3. Relative contribution of FSR versus ISR in the annihilation process e + e~ — ■> 
7r + 7r~7 at ^fs = 1.02GeV for different photon and pion angular cuts. E 1 > lOMeV. 

the main issues of the radiative return method is the suppression of this kind 
of events by choosing suitable kinematical cuts. From EVA studies this is 
achieved by selecting events with the tagged photons close to the beam axis 
and well separated from the hadrons which reduces FSR to a reasonable 
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Fig. 4. NLO corrections to initial state radiation in the annihilation process e + e — > 
7 + hadrons. 

limit. This is illustrated in figure 3. Furthermore, the suppression of FSR 
overcomes the problem of its model dependence which should be taken into 
account in a completely inclusive measurement [20]. 

Preliminary experimental results using the radiative return method have 
been presented recently by the KLOE collaboration [21]. Large event rates 
were also observed at BaBar [22]. 

The theoretical description of the radiative events under consideration 
to a precision better than 1% requires a precise control of higher order 
radiative corrections. The next sections are devoted to the calculation of 
the NLO corrections to ISR and its inclusion in a new improved Monte 
Carlo generator [24, 25]. 



where the virtual photon decays into a hadronic final state, J*(Q) — ► hadrons, 
and the real one is emitted from the initial electron or positron, receives con- 
tributions from one-loop corrections and from the emission of a second real 
photon (see figure 4 for a schematic representation). 

The calculation of the different contributions proceeds as follows: first 
the one-loop diagrams are reduced, using the standard Passarino-Veltman 
procedure [23] , to the calculation of a few one- loop scalar integrals which are 
calculated in D = 4 — 2e dimensions in order to regularize their divergences. 
The interference with the lowest order Feynman diagrams is calculated and 
the ultraviolet (UV) divergences are renormalized using the on-shell mass 
scheme. The remaining infrared (IR) divergences are cancelled by adding 
the soft contribution of the second photon calculated analytically by inte- 



4. NLO corrections to ISR 



At NLO, the e + e annihilation process 




(6) 
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gration in phase space up to a soft energy cutoff E y < w^fs far below the 
center of mass energy y/s. 

In order to facilitate the extension of the Monte Carlo simulation to 
different hadronic exclusive channels the differential rate is cast into the 
product of a leptonic and a hadronic tensor and the corresponding factorized 
phase space 

1 dQ 2 
do = —L^H» u d<S> 2 (pi,P2]Q,k 1 )d<S> n ( y Q- 1 q 1 ,-,q n )— , (7) 

where d$ n (Q;qi,-,q n ) denotes the hadronic n-body phase-space including 
all statistical factors and Q 2 is the invariant mass of the hadronic system. 
The physics of the hadronic system, whose description is model dependent, 
enters only through the hadronic tensor H^ u . 

The leptonic tensor which describes the next-to-leading order virtual 
and soft QED corrections to initial state radiation in e + e~ annihilation has 
the following general form: 



L 



2 



virt+soft Q4 yi y2 



P ^P 2 +yqp 1 p\p 2 - p^p ± 

+ Ol2 r «tt a_i 



(8) 



where yi = 2k\ ■ Pi/s with p\ {p 2 ) the four momentum of the positron 
(electron). The scalar coefficients and a_i allow the following expansion 

<Hj = 4f + - a if . o_i = - a a J , (9) 

where give the LO contribution. The NLO coefficients and a^ 1 ] 
were calculated in [24] for the case where the observed photon is far from 
the collinear region. The extension of these results to the forward and 
backward regions is commented in Section 6. 

Finally, the contribution from radiation of two real hard photons, E~ f > 
Wy/s, is calculated numerically using the helicity amplitude method [26]. 
The sum of the virtual plus soft corrections to one single photon events and 
the hard contribution of two photon emission gives the final result at NLO 
which is independent of the soft photon cutoff w. 



5. Monte Carlo simulation 



A Monte Carlo generator has been built which simulates the production 
of two charged pions together with one or two hard photons and includes 
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virtual and soft photon corrections to the emission of one single real photon. 
It supersedes the previous versions of the EVA [17, 24] Monte Carlo. Again 
the program exhibits a modular structure which preserves the factoriza- 
tion of the initial state QED corrections. The simulation of other exclusive 
hadronic channels can, therefore, be easily included with the simple re- 
placement of the current (s) of the existing modes, and the corresponding 
multi-particle hadronic phase space. The simulation of the four pion chan- 
nel [19] will be incorporated soon as well as other multi-hadron final states. 
Our results will be presented in [25]. 



6. Tagged or untagged photons 

Prom the experimental point of view it would be much easier to perform 
the analysis with no lower photon angular boundary. The cross section in 
the forward and backward regions grows very fast and therefore a small 
deviation in the determination of the photon angle could introduce a large 
error. It turns out that further advantages appear if only an upper cut 
on the photon angle is imposed. The cross section is thus larger. For 
instance, the cross section for radiative events with # 7 < 21° and 55° < 6 W < 
125° is 4 times bigger than the corresponding cross section for 5° < # 7 < 
21°. Furthermore, since FSR is isotropic with respect to the beam axis its 
contribution does not change much and, as a result, its relative importance 
with respect to ISR is much smaller (see lower line in figure 3). All together 
may provide a better control of the systematics of the measurement and, 
therefore, a better determination of the hadronic cross section. 

It should be point out that in the forward and backward regions some 
corrections of order m 2 /s, where m e is the electron mass, has to be taken 
into account even though m 2 /s is a small quantity. At LO, the full leptonic 
tensor is given by: 



Iff 



+ 



8m 2 



2mV(l 



,2\2 



4q 2 



2/12/2 



A 4 v 

ViPi 



2/12/2/ ■- 

P^+PlPt 



8m 2\ ^J.^ 



+ 



2q 2 +y\+vl 



2/12/2 

8m 2 Aq 2 \ p^p v 2 
yj 2/12/2/ s 



(10) 



with yi = 2k\ -pi/s, m 2 = m 2 /s and q 2 = Q 2 /s. Expressing the bilinear 
products yi by the photon emission angle in the center of mass frame 



1 



2/1,2 



-(l=F/?cos0) , = \A - 4m 2 . 
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The electron mass enters twofold. Through j/j, i = 1,2, but also quadrat- 
ically in terms of the type m 2 /yf. In the former case, the electron mass 
regulates the collinear divergences of the cross section. The second kind of 
contributions gives only finite corrections which, nevertheless, can be size- 
able. In both cases the electron mass plays a relevant role only when the 
photon is emitted at very small angles. Otherwise, the limit m 2 /s — > can 
be taken safely. 

To illustrate this point the differential cross section for the process 
e + e~ — > 7T + 7r _ 7 at -y/s = 1.02 GeV is presented in figure 5. The pion 
and photon angles are integrated without applying any angular cut. The 
electron mass is kept in the scalar bilinears yi but the quadratic terms are 
considered or not to show their effect. Both results are compared in the 
small insert of this figure with the known corresponding analytical differen- 
tial distribution 

~ da 4a 3 A s 2 + Q 4 ( s \] , . 
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Fig. 6. Electron mass corrections for different photon angular cuts. 



with the substitution R(Q 2 ) -> (1 -Aml/Q 2 ) 3 / 2 \ F 27T {Q 2 ) \ 2 /4 for the two 
pion exclusive channel, being Fi-k the pion form factor. In (11) the expansion 
in the electron mass has been performed after integration. Figure 5 shows 
thus the good performance of the Monte Carlo integration also at very 
small photon angles. The analytical result agrees with the Monte Carlo 
integration when the full electron mass dependence is considered in Eq.(10). 
If the quadratic mass terms are not taken into account a deviation of up 
to 8% is found. In contrast, at higher energies these corrections are quite 
suppressed, less than 1 per mil at = 10 GeV. 

Figure 6 shows the effect of the electron mass corrections for different 
angular cuts at \fs = 1.02 GeV. Up to 10% deviation with respect to the 
full result is found which nevertheless decreases quite steeply if the collinear 
region is excluded: to less than 5 per mil for tagged photon events at 1 
degree minimal angle. 

These electron mass finite corrections, which are relevant only in the 
very forward and backward regions, below 1 degree at ->/s = 1.02 GeV, have 
been included in the new version of the EVA Monte Carlo. At present, only 
at LO. At NLO, the corresponding electron mass corrections to the virtual 
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and soft contributions in one photon events are included in the definition 
of the j/j's while the quadratic terms will be included soon. For two photon 
events the full electron mass dependence is already considered. 



The radiative return method is a competitive method to measure the 
hadronic cross section at high luminosity e + e~ colliders (4> and -B-factories) 
giving access to a wide range of energies, from threshold to the center of 
mass energy of the collider. A Monte Carlo generator has been built which 
simulates the production of two charged pions together with one or two hard 
photons and includes virtual and soft photon corrections to the emission of 
one single real photon. Its modular structure is such that the simulation of 
other exclusive hadronic channels can be easily included. The description 
of events with hard photons radiated at very small angle has been also 
investigated. 
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